Quantum manipulation via atomic-scale magnetoelectric effects 
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Magnetoelectric effects at the atomic scale are demonstrated to afford unique functionality. This 
is shown explicitly for a quantum corral defined by a wall of magnetic atoms deposited on a metal 
surface where spin-orbit coupling is observable. We show these magnetoelectric effects allow one to 
control the properties of systems placed inside the corral as well as their electronic signatures; they 
provide alternative tools for probing electronic properties at the atomic scale. 



It has long been appreciated that interesting and 
unique properties arise from the couphng between the 
charge and magnetic degrees of freedom in materials. 
Advances in our understanding of the physics regulat- 
ing these properties have given rise to systems with a 
number of practical applications. Prominent examples 
include multiferroic materials, which exhibit simultane- 
ous and cooperative ferroelectric and magnetic orderingji 
as well as giant^ and colossal'^ magnetoresistance materi- 
als, which exhibit substantial changes in their electronic 
transport with the application of small magnetic fields. 
These and other magnetoelectric (ME) effects — cou- 
pling charge and magnetic degrees of freedom — provide 
unique functionalities and hold promise for novel device 
applications.^ 

With the continuing drive to miniaturize electronic de- 
vices, there is interest in better controlling or/and en- 
hancing the functionality of nanoscale systems; there is 
particular interest in novel devices that utilize phenom- 
ena inherent/unique to these nanometer scalesi^ Here, we 
predict unique functionality in a nanoscale device aris- 
ing from the coupling of charge and magnetic degrees of 
freedom in the ultimate miniaturization, namely where 
devices are built atom-by-atom:^-- We demonstrate that 
the interplay of spin-orbit coupling and electronic scat- 
tering results in ME effects which enable exquisite con- 
trol at the atomic scale. This control provides a powerful 
tool for manipulating the response of quantum systems, 
adding desirable functionalities not only for fundamen- 
tal studies, but also for possible future devices built in a 
"bottom up" approach. 

Our device system consists of a quantum corral (QC) 
on a metal surface with spin-orbit coupling (SOC) (e.g. 
Au(lll)^), where the QC's wall is made of magnetic 
atoms. We demonstrate the possibility of controlling the 
electronic properties of the QC by changing the magneti- 
zation of the atoms forming the QC's wall; we show that 
these ME effects allow one to control the properties of 
systems placed inside the QC as well as their electronic 
signatures. This control provides powerful alternative 
tools for probing and manipulating electronic properties 
at the atomic scale. 

The Hamiltonian for the system has the form 
H=Hqc+Hi, where Hqq describes the QC, and Hi de- 
scribes a system we place inside the QC, whose prop- 
erties will be controlled and/or probed (see below). We 



describe the QC by the Hamiltonian Hqc—Hq+V, where 
Ho describes a two-dimensional electron gas (2DEG) 
with SOC, and y is a scattering potential describing the 
QC's wall. The Hamiltonian for the 2DEG is 



^" " 2^^^ + A z • (p X a) 



(1) 



where p is the momentum operator of the 2DEG, {cr^} 
are the Pauli matrices, m* is the electron's band mass, 
and A parameterizes the SOC. As described above, we 
are interested in the case where the QC's wall is made of 
magnetic atoms — being interested in the system's low- 
energy properties, we treat the atoms as a collection of 
s-wave scatterersjS we take 
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where r is the position operator of the 2DEG, Vq de- 
scribes the potential scattering, J is the exchange cou- 
pling between the (magnetic) atoms and the 2DEG, Ti 
is the spin operator of the i*'' atom, and the atoms are 
located at the positions {ri}.- The form(s) of Hi will be 
specified subsequently. 

In this work, we will be interested in the case where 
the atoms of the QC's wall are ferromagnetically (FM) 
ordered — we assume the atoms' moments are sufficiently 
large and treat them as classical variables: (J/2) (Ti)=M. 
The physical quantity of interest is the electronic local 
density of states (LDOS) in the QC, A(r, w). [The differ- 
ential conductance measured in scanning tunneling mi- 
croscopy is proportional to A(r, aj).l°] This is obtained 
from the QC's retarded Green's function (GF) G(r, r'; cu) 
via 



A{r,uj) = -- Im{Tr[G(r,r;a;)]} 



(3) 



when 
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G(r, r'; u) = Go(r, r'; u)+Goir, Ro; uj)Tiu)Go{Ro, r'; lo) . 

(4) 
In Eq. m the T-matrix T{uj) describes the influence of 
Hi. Furthermore, Go(r,r';a;) is the bare GF of the QC, 
i.e. the GF in the absence of Hi^ 



In what follows, we choose (7ryOo)Vb=0.3 and 
(7r/cio)|M|=0.5 {po—m*/2n)-r- we consider physically rea- 
sonable values of the parameters for the 2DEG>i^ a Fermi 
energy Ep=0.5eV, ■m*=0.26me {me is the bare electron 
mass), and A=4xlO~^^eV-m. The results we show are 
for an elliptical QC with 40 atoms, similar to what has 
been realized experimentally^^ {x/a)'^+{y/b)'^—R'^ with 
i?=57.22A, a/6=1.5, and (±c,0)=(±Va2 - 52, 0) being 
the ellipse's foci; we will also comment about results ob- 
tained for a circular QC. It should be stressed the results 
we report are robust — changes in the relative size and 
phase of the ratio M/Vq have only quantitative effects 
on the results, leaving our overall discussion and conclu- 
sions unaffected; furthermore, the corral's geometry can, 
in fact, be tuned to enhance/optimize the ME effects (de- 
pending on the parameters). 

We begin by discussing the electronic properties of the 
QC (with Hi=0). Fig. [T] shows a spatial scan of the 
LDOS at Ep for different directions of M. Notice how 
the LDOS changes as one changes the direction of M — 
the magnetization of the QC's wall and the SOC give rise 
to strong ME effects at the atomic scale. These effects 
are due to the breaking of SU(2) spin-rotation invariance 
by the SOC. Furthermore, the differences in the LDOS 
for M=|M|x and M=|M|?; in the elliptical QC are due to 
the breaking of rotational invariance — for a circular QC, 
the DOS for M=|M|x and M=|M|y are identical. Fig.[l] 
shows the energy dependence of the DOS at particular 
points in the QC — changing the direction of M changes 
the energy dependence of the DOS; indeed, by carefully 
choosing the position in the QC, the changes can be quite 
pronounced [see Fig. [2Jb)]. For reference, the DOS with 
|M|=0 is also shown, as well as the DOS with A=0 (in 
the inset). 

The ME effects exhibited by this system could find 
utility in a variety of applications that exploit the spatial 
and/or energy dependence of the QC's LDOS. To illus- 
trate the utility of the spatial dependence of the LDOS, 
we place a spin-1/2 magnetic impurity at a position Rq 
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FIG. 2: Energy dependence of the LDOS for different M: (a) 
r=(0,0) (b) r=(-c,0). Inset: LDOS at (-c, 0) for |M| = 
with and without SOC. 



inside the QC; wc investigate how its low-energy prop- 
erties depend on its position. We describe the magnetic 
impurity by the Hamiltonian 



Hi = JKr-S(Ro) 



(5) 



where r is the impurity's spin operator, S(Ro) is the 
2DEG's spin operator at the position Rq, and Jk is the 
exchange coupling between the impurity spin and the 
2DEG. In what follows, we take Jk>0-^ 

Interestingly, this seemingly simple system exhibits 
nontrivial behavior in the infrared, due to quantum fluc- 
tuations — a strongly correlated state arises, where a 
cloud of conduction electrons forms a singlet with the 
impurity, U' This strongly correlated state manifests itself 
via a resonance at (or near) the Fermi energy, referred to 
as the Kondo resonance (KR) . This can be seen in Fig. [3l 
where the imaginary part of the impurity's T-matrix is 
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FIG. 1: Spatial scan of the LDOS at Ep: (a) M=|M|t/ (b) 
M=|M|x. 
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FIG. 3: Imaginary part of the T-matrix due to a magnetic 
impurity for different M: (a) impurity at the origin Ro = (0, 0) 
(b) impurity at the focus Ro = (— c, 0). 
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FIG. 4: Spatial dependence of the LDOS difference at Ep due 
to a magnetic impurity, SA{r, Ep) (see text): (a) impurity at 
the focus Ro = (c, 0) (b) impurity at the origin Ro = (0,0). 



shown and, in particular, the KR appears. Furthermore, 
the width of this KR represents the dynamically gener- 
ated scale characteristic of this strongly correlated state, 
the Kondo temperature Tk>^ 

The ME effects allow one to control the Kondo effect 
exhibited by the magnetic impurity in the QC; more gen- 
erally, it allows one to control the magnetic properties of 
an atom placed in the QC. Indeed, Fig. [3] shows that the 
KR can be controlled by changing (the direction of) M. 
Furthermore, we see that the change in the KR depends 
on the position at which the impurity is placed^^ — one 
can control the impurity's properties in a desired way by 
a judicious choice of its position. It should also be noted 
that, besides impacting the properties of the magnetic 
impurity, the QC is also impacted by the magnetic im- 
purity; the impurity's influence on the QC depends on its 
position. This can be seen in Fig. |4l where a spatial scan 
of the quantity 5A{Y,Ep)=A{r,Ep)x—A{Y,Ep)y — the 
difference in the LDOS (at Ep) between M=|M|a; and 
M=|M|y — is shown for different positions of the mag- 
netic impurity. [Results obtained for a circular QC are 
qualitatively similar to those obtained for an elliptical 
QC] 

As we saw above, the ME effects allow manipulation 
of both the spatial and energy dependence of the QC's 
LDOS. As we now demonstrate, this change in the energy 
dependence of the DOS can be used for signal filtering. 
To this end, we place a molecule in the QC at the focus 
(c, 0), such that a vibrational mode (VM) of this molecule 
couples to the 2DEG of the QC; we investigate the image 
or "mirage" of the VM at the other focus (— c, 0). We 
describe the VM by-^^ 



Hi ^ -p^ + -n'^x'^ +g X n(Ro) , 



(6) 



where x (p) is the position (momentum) operator of 
the VM, J7 is its characteristic frequency, n(Ro) is the 
2DEG's density operator at the position Ro=(c, 0), and g 
describes the coupling between the VM and the 2DEG. In 
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5: Maps of \duj5A{r;uj)\ due to the vibrational mode for 

M\x and M=|Mjy: (a) and (b) at r=(c, 0)=Ro; (c) and 

(— c, 0), demonstrating the filtering of the mirage 



what follows, we assume the coupling between the 2DEG 
and the VM to be weak.- 

The QC enables signals to be transmitted between 
foci;^"^ it also provides a "cloak of invisibility", similar 
to what has been achieved with electromagnetic fields,— 
where objects were made invisible within a certain fre- 
quency band. More specifically, the electronic properties 
of the QC govern the energy regime in which a signal 
transmitted from one focus can be observed at the other 
focus; in particular, signals within a certain frequency 
range can be hidden from observation.— This cloaking 
can be seen in Fig. [SJ where a density plot of the en- 
ergy derivative of the LDOS \dujA{r; w)| is shown — even 
when the VM is visible at the focus in which it is sitting 
(Fig. [5ja) and (b)), the QC cloaks it from observation 
at the other focus for fl within a certain range (Fig. [SJc) 
and (d)).^^ In this system, the ME effects allow the range 
over which cloaking occurs to be controlled by changing 
the orientation of M — compare the results for M=|M|a; 
and M=|M|y. Said in another way, the ME effects allow 
one to filter the signal transmitted from one focus to the 
other. 

Up to now, we have considered the system's electronic 
properties, assuming the QC's wall to be FM ordered. 
We have also considered the magnetic properties of the 
wall and, in particular, the wall's magnetic ordering 
tendencies;^ as the QC's wall is one-dimensional, fluc- 
tuations will suppress ordering, and an external field h 
is necessary to stabilize the order. Fig. [6] shows results 
for the QC's magnetization for various values of h. We 
see that the moments are disordered at zero field and, 
hence, the magnetization is zero; a nonzero magnetiza- 
tion is obtained as h is increased. We have found that a 
FM aligned wall is obtained from readily accessible mag- 
netic fields — |h|~0.27r(po'/)^£'F; for reasonable values of 
the parameters, this gives |h|=C'(lmeV). [It is worth not- 
ing such values are considerably lower than Kondo tem- 
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FIG. 6: Evolution of the magnetization with the external 
magnetic field h: (a) |h|=0 (b) lhl=0.04£o (c) |h|=0.12£o 
(d) |hi=0.24£;o where Eo=Tr{poJfEF. 



peratures that have been observed from atoms/molecules 
on surfaces;i2,] Furthermore, we have found the geometry 
can, in fact, be optimized, so that ferromagnetic ordering 
occurs at extremely small values of h. 

This work demonstrates proof of principle of the func- 



tionality afforded by ME effects at the atomic scale; 
indeed, we were able to control the properties of sys- 
tems placed inside the QC as well as their electronic sig- 
nals/signatures. With a FM aligned wall, the ME ef- 
fects allowed us to control the magnetic properties of 
atoms placed inside the QC, as well as to filter trans- 
mitted signals; different magnetization patterns for the 
wallj2£ as well as different QC geometries, could pro- 
vide further flexibility and control. By placing several 
atoms/molecules inside the QC, one could engineer de- 
vices where the ME effects allow the entanglement^^ be- 
tween atoms/molecules to be manipulated. It should also 
be mentioned that the ME effects allow manipulation of 
the properties of the QC's wall — the QC's wall itself 
provides a unique magnetic system with interesting prop- 
erties, which could also afford means of transmitting and 
manipulating information.— i^ 
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Quantum manipulation via atomic-scale magnetoelectric effects: 
Supplementary Information 

I. THE SYSTEM AND HAMILTONIAN 

We consider a quantum corral made of magnetic atoms on a metallic surface with spin-orbit coupling (SOC). The 
Hamiltonian is Hqc^Hq+V, where Hq describes the two-dimensional electron gas (2DEG) of the surface, and V 
describes the coupling of the 2DEG to the magnetic atoms. In second quantized form, the Hamiltonian for the 2DEG 



Ho= j dvi^\v) 



^p2 + A£.(pxa) 



i^(j) , (7) 



where V'^(r) is a two-component field operator for the 2DEG tjj\v) — (■0l(r),-0l(r)), and {cr^} are the Pauli matrices; 
the coupling of the 2DEG to the magnetic atoms is 

V^Y. V^^(r') (^0 + ^r, • a\ ^{v,) , (8) 

where r^ is the spin operator for the magnetic moment of the i"^ atom, Vq describes the potential scattering, and J 
is the exchange coupling between the 2DEG and the magnetic atoms. [As before, {cr^} are the Pauli matrices.] 

As we are interested in the case where the wall is ferromagnetically ordered, we treat the magnetic moments of the 
atoms as classical variables; then (J/2)(ri)— ^M. The quantity entering in the calculations is (7r/Oo)|M|. To estimate 
this quantity, we consider, for concreteness, a spin-5/2 moment (which is relevant to e.g. Mn atoms); we consider 
the physically reasonable value poJ=0.2 — we obtain (7rpo)|M|~0.785. Motivated by this value, in our calculations 
we used (7rpo)|M|=0.5. As noted in the text, however, our results are robust, as the corral's geometry can be tuned 
to enhance/optimize the magnetoelectric effects. 

II. SCATTERING FORMALISM 

The QC's GF can be written as 

G(r, r'; uj) = Go(r, r'; uj) + Go(r, Rq; uj)f{uj)Go(Ro, r'; uj) . (9) 

In Eq. [21 the T-matrix Ticu) describes the influence of Hi. Furthermore, Go(r, r';w) is the bare GF of the QC — it 
is the GF in the absence of Hi ; it is determined by the Dyson equation 

Go(r,r';^) = Goo(r,r';c^) (10) 

+ Yl ^oo(r, i-i; w) {Vo I + M-a) Go(r,, r'; oj) 

i 

with Goo(r,r';w) being the free-particle GF i.e. the GF in the absence of the QC's wall (and also Hi). Goo(r, r';^) 
is given byii^ (for rT^r') 

^uuvi , i ) ^; — "-^0 (-^5 ^) I + Gi {R; w) ( ^ ^^ 
where 



Goo(r,r';..) = GS]°(i?;c.)/ + Gr(i?;-)K3^;(^^) -^.^^v-^^y ) (^^^ 



.m 



Xm\ ^^ .„,, , ,, /^ \m\ 



G°°{R;lo) - -ijUl+'-^\Ho[Rik + \m)]+ll-'-j^\Ho[Rik-Xm)]^ (12a) 

G°0(i?;c^) = -!^ Ul + ^\ Hi[R{k + \m)] - fl - ^^ Hi[R{k ~ Xm)]\ (12b) 

with Ho{x)—Jo{x)+iNo{x), Hi{x)^Ji{x)+iNi{x), and exp(i6') = [(a; — x') +i{y — y')]/\r — r'\. In the above equations, 
Jo (a;) and Ji{x) {No{x) and A^i(a;)) are the Bessel function (Neumann function) of order-zero and order- 1, respectively^ 
Furthermore, R= |r — r'| and k is such that fc^/2m — to + {Xm)'^/2m. 



III. MAGNETIC IMPURITY 

We consider a spin-1/2 magnetic impurity placed in the QC; the Hamiltonian is 

Hi=Jkt- ^^(Ro) {a/2) V^(Ro) (13) 

where r is the impurity's spin operator. [Jk>Q-] To proceed, we treat t with a fermion representation — we write 
T~{l/2)f^Wf where /^ is the two-component spinor f^—{fl, /|); the /-fermions satisfy the constraint f^ f=l. Then, 
Eq. [13] can be written as 

i^i--^(V'URo)/)(/VXRo)) . (14) 

As we are interested in the infrared fixed point of Eq. [13] (or, equivalently, Eq. [T4)) . we employ mean- field theory^ 
— the infrared properties are determined by the effective Hamiltonian 

i?cff = xpf + x V^^(Ro)/ + X /^V^(Ro) , (15) 

where x and A are (constants) determined self-consistently via 

-^X = -(V'^(Ro)/ + /V(Ro)) , (/V) = l • (16) 

Then, the T- matrix in Eq. [5]is proportional to the /-fermions' retarded GF:** T{^) — X^ G'{uj), where 

G/(..) = -ze(i)({/(i),/t})(^). (17) 

IV. VIBRATIONAL MODE 
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FIG. 7: QC's GF, taking into account the VM. The single solid (wavy) line denotes the QC's (VM's) bare GF. 

Assuming the coupling between the VM and the 2DEG to be weak (i.e. g is small), we take the VM into account 
perturbatively. To lowest nontrivial order, the QC's GF (taking into account the VM) is given by the diagram in 
Fig. [71 Performing the Matsubara sum and analytically continuing to real frequencies, the T-matrix (in Eq. ^ is 
given by 

2 [d^ . ^ ,{[l + n{U)-f{v)] ^ [n{n) + f{v)] 



as^ 



T{u:)=g' -A,{v)l ^ ^ ''\'' + ' ' VoV A ' (^^^ 

where Ao{i') is the corral's spectral function at the site Rq — ^o(i^) is such that the Matsubara GF can be written 

diy Ao{h') 
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V. MAGNETIC PROPERTIES OF THE QUANTUM CORRAL'S WALL 

As we are interested in the magnetic properties of the QC's wall, we return to Eqs.[7|and[S| — we integrate out the 
2DEG and obtain an effective spin Hamiltonian;^ the leading interaction generated-^ is given by^i^ 

i^spin = ^J-^ / d^ /(^) Ini{ Tr[(ri ■a)Go{ri,rj;u){Tj • ct) G'o(r-,-, r^; w)] } , (19) 



where /(w) is the Fermi function, and Go{ri,rj;uj) is the retarded GF of the 2DEG in the absence of the exchange 
couphng J. To proceed efficiently, we employ a commonly used approximation, namely approximating Go{ri,rj;uj) 
in Eq. [T9]by Goo{ri,rj;uj)~^ We have checked that along the corral's wall, Go(ri,rj;a;) is similar to Goo(ri, r^; w). 
[Of course this approximation fails inside the corral.] We obtain 

-ffspin = ^ Ki] Ti ■ Tj + Dij Clij ■ {ji X Tj) + J°- (Qy • Ti) (Qij • T j) , (20) 

i<j 

where Qy=zx Rij with i?ij=(Ri — Rj)/|Ri — Rj| (z is the unit vector perpendicular to the 2DEG), and the couplings 
are given by 



K.. 



Dy = 



27r 
j2 rEF 



^Llm{[G™(|r,-r,|;c.)]'+[G?0(|r,-r,|;..)]'} 

duj Re{G°"(|r,; - v,\-uj)Gf{\v, - v,\-lo)] , 

jO ^ -li rLl^{[Gf{\v^~v,\-u)n , 
■"■ Jo '^ -' 



where Ep is the Fermi energy. [GQ^dr^ — Vj\\uj) and G^{\Yi — rj|;aj) are given in Eqs. 112a] and I12bl ] 

The magnetic properties of the QC's wall were determined by Monte Carlo simulations of Eq. [501 — the magnetic 
moments were treated as classical three-dimensional vectors of unit length; the minimum energy state of Eq. [20] was 
obtained via a simulated annealing procedure. 
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